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Abstract. We theoretically analyze the Bragg spectroscopic interferometer of two 
spatially separated atomic Bose-Einstein condensates that was experimentally realized 
by Saba et al. [Science 2005 307 1945] by continuously monitoring the relative 
phase evolution. Even though the atoms in the light-stimulated Bragg scattering 
interact with intense coherent laser beams, we show that the phase is created 
by quantum measurement-induced back-action on the homodyne photo-current of 
the lasers, opening possibilities for quantum-enhanced interferometric schemes. We 
identify two regimes of phase evolution: a running phase regime which was observed in 
the experiment of Saba et al, that is sensitive to an energy offset and suitable for an 
interferometer, and a trapped phase regime, that can be insensitive to applied forces 
and detrimental to interferometric applications. 



PACS numbers: 03.75.Dg,03.75.Gg,37.25.+k,67.85.Jk 
1. Introduction 

Bragg spectroscopy has become an established spectroscopic tool in ultracold atom 
experiments [TJ |2j [3j HJ El E] . In typical set-ups an intersecting pair of low-intensity pulsed 
laser beams is used to excite atoms to higher momentum states. The momentum kick 
experienced by the atoms corresponds to the recoil of a photon upon light-stimulated 
scattering between the two laser beams. As the spontaneous scattering for off-resonant 
lasers is negligible and the photons are only exchanged between the directed coherent 
laser beams, the momentum transfer of the atoms can be measured for specific values of 
the energy and the momentum. In particular, in a spectroscopic analysis of the many- 
particle properties of ultracold atoms it is sufficient in the scattering process to describe 
the light beams classically. 

In the experiments by Saba et al. [7] the relative phase coherence between two 
Bose-Einstein condensates (BECs) was measured by Bragg scattering atoms between 
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two condensate fragments. Previous Bragg spectroscopy experiments based on time of 
flight had concentrated on directly detecting the atoms that were transferred to higher 
momentum states by the laser beams. In the experiment by Saba et al. [7J, however, 
the strength of the Bragg scattering was measured by monitoring the variations of the 
light intensity in the laser beams by homodyne detection. Due to the correspondence 
between the light-stimulated scattering of photons between the laser beams and the 
atoms scattered between two momentum states, the intensity fluctuations are directly 
proportional to the number of atoms scattered between the condensate fragments. 




(W,ki) (w 2 ,k 2 ) 



Figure 1. Our model of the Bragg interferometric measurement of the relative 
phase between two distant BECs. The two BECs are described by the macroscopic 
wavefunctions </>&(r) and </> c (r), and are illuminated by two coherent laser beams. Bragg 
scattering imparts momentum to atoms from the left condensate, transferring them 
to the state described by 0fc(r). After an appropriate time the outcoupled atoms will 
overlap with the right condensate, and the Bragg beams will drive Rabi oscillations 
between the two atomic clouds. This establishes an optical weak link between the 
two BECs, and continuous monitoring of the intensity fluctuations in the laser beams 
measures the phase coherence between the BECS. 

Saba et al. [7] measured the light intensity variations of the Bragg beams which 
revealed relative phase coherence between the condensates even when the BECs were 
independently produced and possessed no a priori phase information. By analyzing 
theoretically a continuous atom detection process it has been previously shown that 
the back-action of quantum measurement of the atomic correlations [8j [9j [TOj [TTJ [121 
IT3l HU [15], and analogous photon correlations [16j [17], can establish a relative phase 
between two BECs even when they have 'never seen each other' before. It has also been 
suggested that phase-coherent states of condensates may naturally emerge as robust 
state descriptions due to dissipative interaction with the environment [T8| [15] . With 
regards to the Bragg spectroscopic interferometer of [7J, the question we ask is: how is 
the phase coherence between the two BECs created, given the condensates interact with 
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coherent laser beams that can usually be described classically? 

Here we analyze a model of the experimental detection scheme [7J, illustrated in 
figure [U and show that the phase coherence can be built up by continuously monitoring 
the photo-current obtained by a homodyne measurement that describes the intensity 
fluctuations of the laser beams. We show a rapid establishment of a well-defined relative 
phase between two independently produced BECs. We identify two distinct regimes of 
subsequent phase evolution: a running phase and a trapped phase regime. In the 
running phase regime the relative phase grows linearly proportional to the energy offset 
between the two condensate wells and could be suitable for a weak force detection in 
interferometric applications [7]. In the trapped phase regime, in the case of a very weak 
energy offset, the measurement process drives the system close to a dark state where 
a destructive interference between different scattering paths suppresses the intensity 
fluctuations of the lasers. In the trapped phase regime the effect of the energy offset 
on the phase evolution is suppressed, potentially to the detriment of interferometric 
applications. 

Our analysis demonstrates how Bragg spectroscopy can be sensitive to subtle 
quantum features of ultracold atom systems. Quantum measurement-induced back- 
action of the photo-current detection on the relative phase coherence of the BECs 
represents a spatially non-local entanglement of the laser beams and the relative 
many-particle state of the atoms. Indeed, the location of the photo-current detection 
can be far separated from the interaction region of the coherent laser beams and 
the atoms. Moreover, one Bragg pulse can be used to entangle the two spatially 
isolated BECs. A second pulse may then be employed in optical readout of the 
subsequent evolution dynamics of the measurement-established relative phase coherence 
between the condensates. An energy offset between the two condensate wells between 
the subsequent pulses would result in a detectable phase shift providing potential 
interferometric applications [7]. Here the phase is determined by a continuous quantum 
measurement process opening possibilities for quantum feedback and control methods, 
e.g., in generation of sub-shot-noise phase-squeezed states [201 EJJ- Such states may be 
useful in quantum-enhanced metrology in the realization of a high-precision quantum 
interferometer overcoming the standard quantum limit of classical interferometers 
[221 1231 EU [25]. Probe field response was also recently measured in Bragg spectroscopy 
of condensate excitations in a heterodyne-based detection system, which was able to 
reach the shot-noise limit [26J. Previous theoretical studies of the effects of continuous 
monitoring on light scattered from BECs have considered photon counting [T6| ITT], e.g., 
in preparation of macroscopic superposition states [17], and dispersive phase-contrast 
imaging [27J [28], e.g., in suppression of heating [28] . 

The article is organized as follows. In section [2] we give a short review of the 
experimental setup of [7] and the relevant results for this work. We then introduce 
our basic theoretical model. In section [3] we derive a stochastic differential equation 
which describes the evolution of the system under the continuous measurement of the 
scattered light intensity. In section 0] we present our numerical results with a physical 
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interpretation. Finally, some concluding remarks are made in section [SJ 

2. Model and effective Hamiltonian 

An interferometric scheme between two spatially isolated BECs was experimentally 
realized in [7] without the need for splitting or recombining the two condensate atom 
clouds. The method was based on stimulated light scattering of a small fraction of the 
atoms, only weakly perturbing the condensates and therefore representing an almost 
nondestructive measurement. Two isolated BECs were prepared in the sites of an 
unbalanced double-well potential, and illuminated by the same pair of Bragg beams. 
These beams outcoupled atoms from each well, and interference between such atoms 
provided a coupling between the BECs. When outcoupled atoms from one condensate 
spatially overlapped the second, measurement of the Bragg beam intensity was shown to 
be sensitive to the relative phase $ between the condensates. In addition, the potential 
offset between the two wells gave rise to a difference in energies 8fi, which in turn led 
to a relative phase evolution $(t) = $(0) + 5fit/h. This was observed as oscillations in 
the Bragg beam intensity of frequency u osc = 5/i/h, demonstrating how monitoring the 
Bragg beam intensity directly measured the dynamical evolution of the relative phase 
between the macroscopic wavef unctions. 

In the experiment a single Bragg pulse established a random relative phase between 
the two independently produced BECs. If two successive Bragg pulses were applied to 
the same BEC pair the relative phase measured by the second pulse was correlated with 
that detected by the first pulse, indicating that the interaction of the first Bragg beam 
with the atoms had projected the system into a state with a well-defined relative phase 
between the condensates. 

The key to the method is the weak link established between the BECs by the 
Bragg laser beams that couple out small atomic samples from the condensates [29] . The 
coherently driven population dynamics between the BECs is influenced by the relative 
phase coherence [301 EI] > an d the Bragg scattering may be understood as an interference 
in momentum space [32]. The specific advantage of the Bragg spectroscopic interference 
scheme [7] is the nondestructive nature of the detection process, potentially constituting 
a major advance for interferometric applications since it allows one to probe the evolution 
of the phase coherence in time by a continuous measurement process [331 EH EH] . It has 
also been argued that this setup can be seen as an analog to homodyne detection for 
matter waves [36] . 

In order to analyze the continuous measurement process of Bragg spectroscopy we 
consider the system depicted in figure [T] that is analogous to the experimental set-up 
of [7]. We assume that the two condensates are initially uncorrelated and that there 
is no tunnelling between the two spatial regions. As in the experiment, an offset in 
the trapping potential between the two condensates is accounted for by a difference in 
chemical potential, 5fi. The condensates are illuminated by two Bragg beams, which 
impart momentum, kicking atoms out of the traps. The outcoupled atoms propagate 
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from the left condensate to the right and establish an optical weak link between the two 
macroscopic wavefunctions [29] . 

For simplicity, in the theoretical analysis we use a single mode approximation for 
the condensates and assume that all atoms in the left (right) condensate are in the state 
| b) (|c)). The atoms in the left (right) condensate are then described by the second 
quantized field operators "0l(i") = 0t(r)6 {^r(t) = c ( r )c) which fulfill the usual bosonic 
commutation relations. Here b (c) annihilates an atom in state |6) (|c)) and </>b, c ( r ) obey 
the Gross-Pitaevskii equation [37] . 

Atoms from the BEC in the left well in state \b) are transferred by the Bragg beams 
to the momentum state k = ki — k2, where kj are the wavevectors of the Bragg beams. 
The outcoupled atoms propagate with momentum k towards the right BEC in state |c). 
We take the wavefunction of the outcoupled atoms 0&(r) to be the momentum shifted 
original wavefunction 

</» fc (r) = 6 (r-r L )e ik -( r - r ^, (1) 

where tl (i\r) gives the position of the centre of the left (right) trap. We assume that the 
momentum kick of the atoms is sufficiently strong, so that the essential characteristics 
of the continuous quantum measurement process are not obscured by collisions with the 
remaining trapped atoms, collisions among the outcoupled atoms, and the effect of the 
trapping potential. We assume that enough time has passed such that the outcoupled 
atom cloud from the left condensate completely overlaps the right condensate. We 
hence neglect the time evolution of the outcoupled cloud while flying from the left to 
the right trap. In our model, this evolution leads to an additional phase factor which is 
inconsequential to our findings. We also take the same functional form of the trapping 
potential for the atoms in the left and right condensate such that </>&(r— 1) = c (r) = 0(r) 
where 1 = — 17, is the distance vector between the two potential minima. With these 
assumptions we find for the effective Hamiltonian 

Hefi = Ha + Hal + H E m, (2) 

where 



H, 



^] 6ta + ^ 5 t ii + ^( a » St + i . 6 ) (3) 

describes the Rabi oscillations between the outcoupled atoms and the atoms in the 
right condensate due to the Bragg beams. Here Qj are the Rabi frequencies of the Bragg 
beams, A is the detuning from the excited state |e) which couples the two photon Raman 
transition between \bk) <-> \c) and the operator bu annihilates an outcoupled atom in 
the momentum shifted state \bk) with wavefunction 0&(r). Hamiltonian (j2J) is written 
in the reference frame of the Bragg beams where we assume the two laser frequencies 
to be equal coi ~ U2 = ojl- The term 

Hem = th}, A A a A a A (4) 
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takes into account the electromagnetic (EM) vacuum energy, where we have used the 
standard plane wave decomposition for the EM-field modes. Specifically, the positive 
frequency component of the vacuum electric field amplitude reads 



5E + (r,t) = ^y|^e A a A (t)e ik - r . (5) 

A 

Here A labels a mode of the EM-field at wavevector k,\, polarization 6a _L k^, and 
frequency u\ = c\k\\. The velocity of light is denoted by c, the quantization volume 
by V, the vacuum permittivity is e , and A A = uj\ — ujl- The operator a x annihilates a 
photon in mode A. The total electric field is the sum of the coherent Bragg laser fields 
Ej£ -(r) and the vacuum fields 5E + (r). The coherent part is responsible for the driving 
terms in (j3J), while <5E + (r) provides the coupling of the vacuum modes with the atomic 
dipoles. We consider off-resonant scattering where the scattering rates for sufficiently 
large condensates are proportional to the amplitudes of the macroscopically occupied 
modes due to Bose-enhancement, and we neglect scattering to other motional states of 
the atoms. The coupling between the vacuum modes and the atoms is then given by 

H AL = h^2(a[B x + Bla x ^ , (6) 

A 

where we have introduced the operator 

Bx= ((At)*a 1 + (A x 2 ya 2 ), (7) 

with 

eh = c f e, <7 2 = b\e , (8) 



and, after adiabatic elimination, the excited state annihilation operator can be written 
V A A 

We have defined 



A) = (dj ■ e A ) I dr|*(r)| a e-^-^-' , (10) 

where the factor outside the integral is the coupling strength between the atomic dipoles 
and the electromagnetic field mode A [38] . Here the matrix elements of the dipole 
moment operator d for the transition are denoted by d]~ = (c|d|e),clr = (bk\d\e). 



3. Continuous homodyne measurement 

We consider the condensate and the outcoupled atomic cloud together with the driving 
fields as an open quantum system and eliminate the vacuum EM field modes. The aim 
of our treatment is to compute the evolution of the reduced system under continuous 
measurement of the light intensity of the Bragg beams. The intensity of the beam j is 
given by 

J J - = 2ce (E7(r,t)Et(r,t)). (11) 
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Here the total electric field amplitude of each Bragg beam is given by the sum of the 
coherent driving laser field and the field 5Et(r, t) due to scattering in the direction of 
the beam j 

E+(r,t)=E+ J (r) + < 5E+(r,t). (12) 

Assuming that the amplitude of the scattered field is small compared to the applied 
laser field, the measured intensity is approximately 

Ij ~ 2ce ((Er nj E+ nj ) + (Er/E+) + (*EJE^>) , (13) 

where the last two terms give rise to fluctuations in the intensity incident on the detector 
J- 

We may now solve the intensity fluctuations by calculating the scattered field 
amplitude from the effective system Hamiltonian. From the Heisenberg equation of 
motion for a\ one finds 



a x (t) = a A (0)e" iA ^ - i / dt' e - iA ^ ((A^a^f) + (A^& 2 (t')) , (14) 



with (dx{0)) = 0. Inserting f[T4"|) into §5§ and defining q,, = fc^n — k,, as the change 
of the wavevector of light upon scattering with = ujl/c [39], we then obtain two 
contributions from the vacuum field, one for each beam 

h.2 „ik L r D t- 

5E}(r,t) ~ ^ - n x (n x dj)aj{t) J dr'|0(r')| 2 e iq '" r . (15) 

The spatial integral over the wavefunction 0(r') enforces an approximate momentum 
conservation, so that the photons are dominantly scattered into a cone centred at 
qj = in the direction of the laser beam j. In deriving ([15]) we made the expansion 
|r — r'| = re — n • r', with n being the unit vector that points from the scattering region 
to the detector, r D is the distance between the detector and a representative point at 
the origin of the scattering region. Due to the normalization of the wavefunction we 
finally find for the scattered electric field in the two outgoing beams 

1,2 ik L r D 

5E+(r,t) = x (n x dj^t) . (16) 

The atomic operator associated with the spontaneous emission of a photon into beam j 
in ([TBI) is given by &j. The master equation which describes the evolution of the reduced 
density matrix after elimination of the vacuum field modes then reads [40] 

'p(t) = j: [p, Ha]-J2^ {*l*iP + P*fa ~ 2 *iP*i) ■ 
i=i,2 

Here jj is the rate of spontaneously scattered photons, and is related to the total 
spontaneously scattered light intensity 51 = 2c£ (<5Ej(r, t)5E+(r, t)) via [T7] 

J- j dfotSI^yiafa), (18) 
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where the angular integral is over the scattering cone of beam j. The operators which 
are associated with the light field amplitude of the beam j read 

Cj = + crj), (19) 

where atj is proportional to the amplitude of the coherent laser beam with wavevector kj. 
The intensity is then proportional to (CjCj). The leading contribution comes from the 
coherent intensity oc 7j|<x,'| 2 [corresponding to the first term in ( 1 1 3 p ] and the intensity 
fluctuations are dominated by the terms ■jj(aj(aj) + C.c.) [corresponding to the second 
and the third term in (IT51) ]. Extending the treatment of [H] for our setup one finds that 
the evolution of the system under the continuous monitoring of light intensity can be 
described by the stochastic differential equation 

+ dt)) = ( 1 - jH A dt + J2 \j*fa dt + 2 ^ x J dt + °jV^ dW j\ ) IV>(*)>> (20) 



•V, ~ fa ■ o))- (21) 



where 

Xa = 

2 

Here dWj is a Wiener increment with zero mean (dWj) = and ((dWj) 2 ) = dt, 
which appears as a result of the continuous measurement process. Keeping terms to 
lowest order in the fluctuations one finds an expression for the photocurrent in essence 
equivalent to (TT3]) 

if °\t) = irf + a, (2 7i (X i ) + y/yti®) ■ (22) 

Here 

W) = ^ (23) 



represents Gaussian white noise [H] and arises from the open nature of our quantum 
system. 

4. Numerical Results 

In order to study the effect of the homodyne photo-current measurements on the system, 
we numerically integrate fTSOj) . using the Milstein algorithm [42J. As an initial state in 
the numerical simulations we take a pure number state in each condensate, with no 
well-defined phase between them, and the incident Bragg laser beams are taken to be 
classical coherent states. The relative phase between the condensates as a function of 
time may then be calculated as $(£) = arg((c'6}). We define a measure for the strength 
of the phase coherence between the condensates by the absolute value of the normalized 
phase coherence 

9(t) = -JB^. (24) 

;ctc)(StS> 
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Figure 2. (colour online) Coherence g(t) as a function of time t in units of 
to = 2nH/Sfi. The simulation is done for a total of N = 100 atoms which are initially 
distributed equally between the states \bk) and \c). Parameters were chosen to be 
7i = 72 = 10 5 /£o and -y^if^ = 10 3 /io- The black solid line corresponds to A = IOO71 
and the red dashed line to A = IO71. In the inset we show the evolution of the phase 
as a function of time t in units of to- 



A value of g(t) close to one indicates a high degree of relative phase coherence, while 
condensates with no relative phase information have g(t) ~ 0. 

In figure |2] we plot the time evolution of the coherence and the relative phase for 
two different values of the detuning A. No well-defined relative phase exists at early 
times, the coherence starts at zero and $(£) shows large random fluctuations with time. 
As the continuous measurement proceeds the coherence builds rapidly, leading to a well- 
defined relative phase with a stable value. Once established, we then see two different 
regimes of behaviour at longer times. For large values of the detuning A = IOO71 we 
see a running phase behaviour: once well established with a value which is random for 
each individual run, the phase grows linearly in time with a rate proportional to the 
difference in energies between the condensates ~ 5fit/h. From ( 1221) we note that 
the measured photocurrent from the two Bragg beams is essentially proportional to the 
quadrature (Xj) after subtracting the background current, and the corresponding time 
evolution of (Xj) is shown in figure |3j In the running phase regime, the quadrature 
exhibits well defined oscillations with frequency u osc = Sfi/h, and this corresponds to 
the experimental measurements obtained by Saba et al. [7\. Such oscillations thus 
give an interferometric measurement of the relative phase evolution, sensitive to any 
accumulated phase shift due to an energy offset between the distant condensates. An 



Bragg spectroscopic interferometer and measurement-induced correlations in BECs 10 



interferometer of this type could be used, for example, to detect a weak force applied 
to one of the condensates. 

Choosing a smaller value for the detuning A = IO71 we find a very different long 
time behaviour. Once again the phase fluctuates as coherence is established, although 
this occurs on a much faster timescale. This can be understood from the fact that the 
phase is established as a result of the intensity fluctuations in the laser beams, which 
are enhanced by decreasing the detuning A. Unlike the running phase case, once firmly 
established the phase now locks to an almost constant value near to it. This trapped 
phase state has the two condensates almost entirely out of phase, leading to destructive 
interference in the oscillations between the states \b k ) and |c) and resulting in a state 
analogous to a dark-state. The corresponding quadrature (Xi) therefore exhibits merely 
random fluctuations which would not be suited to an interferometric type experiment. 
Note that while one would expect the amplitude of these random fluctuations to be 
suppressed compared to the coherent oscillations of the running phase regime, and this 
is indeed the case, the two different values of detunings used here do not allow such a 
direct comparison in figure [3j 

The two different regimes of behaviour resemble the AC- Josephson and self-trapping 
behaviours seen in double- well condensates jUJ HU US] , although we emphasize that here 
the coupling occurs due to the non-local measurement process induced by the Bragg 
beams. The trapped phase behaviour is more akin to a dark state, however, due to 
the lack of any nonlinearity in the Hamiltonian which is required for macroscopic self- 
trapping. The different regimes may be understood if we assume that a well defined 
phase and population can be associated with each condensate, and neglect any processes 
other than those included in Ha, leaving a two- mode model similar to that considered 
in [13]. The trapped phase regime then occurs with a stable relative phase difference of 
7r when 

fiifio \z\ , \ 

* = 2S ir7r=J- (25) 

where z = (N k — N c )/(N k + N c ) is the relative population difference, where N k ^ is 
the population in state \b k )(\c)). The trapped phase regime therefore requires either 
5{i/h ~ fi 1 f2 2 /A, or a large population imbalance. This is in agreement with our results, 
where the larger value of detuning has 5ji/h ^> Q1O2/A, and the initial population 
balance is not extreme. The trapped phase condition is then not satisfied and we 
observe a running phase behaviour akin to the AC- Josephson effect. Note, however, 
that the dissipation is vital for establishing the relative phase in the first place, and has 
the effect of shifting the relative phase in the trapped regime away from tt in figure [2j 
The model ( )25|) does not specify the role of spontaneous scattering which determines 
the rate at which the system is driven towards the trapped phase state. 

In the experiment of [7], the system typically contained on the order of 10 6 
atoms, although numbers outcoupled would be only a small fraction of this. Numerical 
simulations in our basis for such large numbers are prohibitively slow, so here we have 
typically used a total atom number of 100. Our results show no significant dependence on 
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Figure 3. (colour online) Mean value of the quadrature {X±) as a function of time t. 
The black solid line corresponds to A = IOO71 and the red dashed line to A = IO71. 
Other parameters are as in figure ([3]) . The inset shows an expanded view of the initial 
behaviour. 



atom number, however, and we expect our results to give a good qualitative comparison 
with the physics exhibited in the experiment. We have so far discussed our results 
in dimensionless units; in order to give a specific example we may take 7^2 to be of 
the order of 2tt x 10MHz, for instance. Our results then correspond to the values 
§H/K ~ 2?r x 630Hz, and y/H^h ~ 2vr x 0.45MHz. 

Parameters used in the experiment of [7] were 5fi/h ~ 2it x 1kHz, A = 2n x 1GHz, 
y/Q^h ~ 2tt x 0.45MHz. This yields a ratio n = 2htt 1 Q 2 /{A5n) ~ 0.4. In order for the 
trapped phase regime to be observed the population imbalance would then be required 
to satisfy z 0.92. In the experiment, the actual population in the momentum shifted 
state (\bk)) that overlapped the second condensate (|c)) was of the order of 2 x 10 4 
atoms during the coupling. The corresponding population imbalance was in excess of 
0.96, and hence did not satisfy condition ( |25l) for the trapped phase behaviour. The 
observed running phase behaviour in the experiment is therefore consistent with our 
model. 



5. Concluding remarks 

Bragg spectroscopy was used in [7j to measure the relative phase between two initially 
uncorrelated BECs. By studying a simplified model containing the essential ingredients 
of the experiment, we have demonstrated how the homodyne measurement process 
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builds up a coherent relative phase between the two condensates. This quantum 
measurement-induced backaction entangles the two macroscopic many-body states even 
though the measurement location can be far removed from the region of interactions. 

Following the establishment of a coherent phase, we have identified two distinct 
behaviours under continual subsequent measurement. With a larger atom-laser detuning 
A, or a large initial energy imbalance 5{i, we reproduce the experimental findings of 
m , with the measured photon flux exhibiting oscillations at a frequency corresponding 
to the energy offset of the separated condensates. In this case, once the coherence 
and a random-valued phase is established it evolves linearly in time $(£) ~ 5fit/h. 
Measurable oscillations in the laser beam intensity mean this state has applications to 
quantum-enhanced interferometry, and the measurement backaction could potentially 
be used further to implement feedback mechanisms [21]. By choosing a smaller atom- 
laser detuning, we found instead that the system stabilized to a trapped phase state 
with the condensate relative phase fixed at almost 7r, while the scattered light intensity 
showed only random fluctuations. A semiclassical model can qualitatively describe the 
difference between these two regimes, with the trapped phase behaviour occurring when 
( |25l) was satisfied. 

When atoms from two initially uncorrelated condensates overlap, we have shown 
that Bragg coupling and continuous homodyne measurement can rapidly establish a 
well-defined relative phase. A closely related experiment [46J has been performed using 
ultra-slow light pulses, in which optical information was coherently transported between 
two spatially separated condensates by a travelling matter wave. An ultra-slow light 
pulse was stopped in the first condensate, creating a dark-state superposition between 
two atomic internal states. Upon stopping the pulse one of the internal states received 
a momentum kick, and outcoupled atoms in this state passed through a second distant 
condensate. By illuminating the second condensate with a coupling laser it was possible 
to revive the initial light pulse, even when the BECs were independently produced. In 
the case where the condensates had been prepared separately, a rapid establishment of a 
coherent phase in a manner similar to that described in this paper explains the recovery 
of the light pulse. 
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